Magnetic frustrations in the face centered cubic antiferromagnet NiS2 
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Neutron scattering experiment on NiS2 single crystal revealed a honeycomb pattern of the intensity 
distribution in reciprocal lattice space (continuous-line structure along the fee zone boundary ) 
providing the first direct evidence for nearly frustrated antiferromagnetism (AF) on the face centered 
cubic (fee) lattice. A small but finite lattice distortion below 30.9 K lifts the degeneracy of the 
magnetic ground state due to the frustration and eventually result in the coexistence of the type I and 
the type II AF LROs, which are mutually incompatible in the fee symmetry at higher temperatures. 

75.25.+Z, 75.40.Gb, 75.50.Ee 



In spite of extensive theoretical and experimental ex- 
plorations [01 of frustrated magnetism, search for such 
new examples is still very active. This class of materials 
does provide many new and unexpected phenomena in- 
cluding the high temperature superconductivity. There- 
fore the research and development for such materials con- 
stantly brings new challenges for theoretical interpreta- 
tion. These materials often show the robust magnetic 
frustration phenomenon. A prototype is a Kagome an- 
tiferromagnet, where the magnetic frustration is caused 
by solely rigid geometrical spin structure with a near- 
est neighbor (NN) antiferromagnetic (AF) interaction on 
a 2-dimensional (2D) triangular lattice Extensive 
studies of this model have revealed that suppression of a 
long range order (LRO) in geometrically frustrated com- 
pounds is directly affected by strong quantum fluctua- 
tions. Another example is the 3D antiferromagnet with 
corner-sharing tetrahedra such as transition metal py- 
rochlores and spinels On the other hand, one can 
introduce a fragile frustration phenomenon which is sub- 
ject to both geometrical spin structure and a subtle bal- 
ance among various magnetic interactions. Such a subtle 
balance is easily broken by structural phase transition 
(SPT) as observed experimentally H,^. 

In order to study the spin dynamics of these systems, 
the momentum and energy dependence of the generalized 
spin susceptibility x"(q, cj) is indispensable. However in 
most cases, the lack of large enough single crystal pro- 
hibits the detailed studies on the momentum dependence, 
and often the parameters in the theoretical model have 
been set rather arbitrarily. It is emphasized here that the 
geometrical structure only is not enough to guarantee the 
frustration and degeneracy. In this paper we report the 



neutron scattering experiment on the NiS2 single crys- 
tal with fee lattice structure. The detailed analysis of 
the scattering intensity provides the first direct evidence 
that the system is very near the phase boundary between 
the two types of antiferromagnetic long range order (AF 
LRO). This quasi-degeneracy is lifted by the SPT at low 
temperatures. 

NiS2 possesses anomalies in its thermodynamic prop- 
erties [^-^ and an unusual coexistence of type I and type 
II antiferromagnetic (AF) LROs which are incom- 
patible from symmetry considerations for the undistortcd 
fee lattice. It shows magnetic phase transitions at low 
temperatures with persistent short range magnetic cor- 
relations up to high temperature in consequence of the 
frustrated magnetic structure of NiS2 [1). In addition, 
the structural instability pertinent to this compound was 
clearly observed both in dilatometric and X-ray 

diffraction [|ll].^ measurements. 

Extensive investigations for a series of transition metal 
dichalcogenides, MX2, during the '70s |]l3| appeared to 
conclude that a simple view of the narrow band con- 
cept applied for the d electrons on M sites is sufficient 
to explain most of the experimental facts. However, the 
recent photoemission experiments indicate that NiS2 is 
not even a simple Mott Hubbard insulator but rather a 
charge transfer insulator (CTI) jl^. Note that such CTIs 
are commonly realized in the strongly correlated electron 
systems of the progenitors of high Tc superconductor and 
colossal magneto-resistance materials. 

Experimental study was done using 9 single crystals 
grown by the chemical vapor transport method, with 
a total volume of 0.9 cm^, which were aligned in the 
(hkk) reciprocal lattice zone using shiny (100) and (111) 
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facets. Neutron scattering measurements were performed 
on TOPAN spectrometer installed at the JRR-3M Reac- 
tor of the Japan Atomic Energy Research Institute. The 
crystal structure was determined to be of high symme- 
try (Pa3) above Tn2- Although a small distortion occurs 
below T/V2, it is ignored in our notation of the crystal 
symmetry presentation and we treat all the neutron scat- 
tering data within the fee unit cell notation. 
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FIG. 1. The temperature dependence of the peak intensi- 
ties of the type I (a) AF Bragg (100) and the type II (b) AF 
Bragg (1/2, 1/2, 1/2) peaks. 

Elastic neutron scattering measurements confirmed the 
existence of [0,^ two types of sharp Bragg reflections: 
the type I LRO with AF wave vector Q/ — (1,0,0) ap- 
pears at Tmi — 39. 6K and the type II reflections with 
Q// = (1/2, 1/2, 1/2) arise abruptly below Tn2 = 30.6 K 
superposed on the type I AF LRO with rather discontin- 
uous phase transition (Fig. 1). Furthermore, evaluation 
of the magnetic entropy estimated from the data of heat 
capacity reveals that the sum of the entropy released at 
both transitions (A5 is 0.51 and 0.60 (J/mole-K) for Tjvi 
and Tn2i respectively) is much smaller than R ln3 (9.13 
J/mole-K). In order to probe spin fluctuations which is 
characteristic of these unusual properties of the phase 
transitions, we have performed inelastic neutron scat- 
tering measurements with momentum transfers in the 
(/lOO) — (Ofcfc) plane of reciprocal space where both the 
type I and the type II AF Bragg reflections are observ- 
able. 



~~~~~~.„Jemperature| 

^^^^ (A) -[^=z 


0.6 K (^) T,^i=; 


9.6 K (C) T~150K T~300K 


LRO and 
spin wave 


type 1 
+type II 


type 1 


no Bragg reflection 


diffuse scattering 


Q-(1,0,0) 
(type 1) 


inelastic 
ti)~12 meV 


inelastic 
ti)~12 meV 


quasi-elastic ^ 
scattenng 


Q~4J1) 
u (2,2,2' 

(type II) 


inelastic 
ti)~12 meV 


auasi-elastic ^ 


zone 
-boundary 


none 


scattering 



TABLE I. Three characteristic temperature regimes of 
magnetic excitations in NiS2. 



Spin fluctuations in NiS2 were found to be dramati- 
cally different in three temperature regions as shown in 
Table 1. (A) In the first regime of the paramagnetic 
phase above T^ri, strong antiferromagnetic spin correla- 
tions persist up to high temperatures. Spin fluctuations 
in this phase are manifested by two components: (i) mag- 
netic diffuse scattering around the type I AF Bragg point 
and (ii) magnetic diffuse scattering component along the 
fee zone-boundary around the type II AF Bragg point. 
Both components are observable even at room tempera- 
ture but grow below 150K, where anomalies in the con- 
ductivity appear. (B) In the second, intermediate regime 
Tn2 < T < T/vi, type I AF LRO arises and concomi- 
tantly spin-wave excitations occur. Two diffusive mag- 
netic scattering components still exist as well. However, 
the diffuse scattering around the type I AF Bragg point 
becomes rather inelastic as described below. (C) In the 
low temperature regime below Tn2 , where both the type 
I and II AF LROs are present with a possible rhombohe- 
dral lattice distortion, the spin waves for each AF LRO 
were observed. On the other hand, the diffuse scattering 
around both the type I and II AF Bragg points becomes 
inelastic and the zone boundary component disappears 

Next, we present the detailed experimental results for 
the intermediate regime B followed by theoretical analy- 
sis. At T = 32 K, we measured the Q-scans for the flxed 
energy transfers of w = 2 meV and 10 meV along the 
[100] direction around the (0, 1 -I- Ci 1 + C) points with 
C = 0,0.125,0.25,0.375,0.5 (Fig. 2a). For energy trans- 
fer w = 2 meV (open circles), a well defined single peak 
was observed for < C < 0.25, whereas it splits into two 
distinct peaks for 0.25 < C < 0.5. Note that the peak for 
C = coincides with the type I AF Bragg point (0, 1, 1) 
and that both (±1/2,3/2,3/2) peaks are at the type II 
AF Bragg points. These peaks are located on the zone- 
boundary line of the fee lattice as shown by broken line 
forming a honeycomb pattern in Fig. 2a. The same set 
of scans with lj = 10 mcV (closed circle) reveal basically 
the same geometric picture, although the signal intensity 
is more concentrated around the type I AF Bragg point. 
(We define this as a type I peak.) The neutron scattering 
signal is purely magnetic, since the scattering intensities 
at momentum transfer (300) follows the square of the 
Ni^+ ion magnetic form factor. 

We observed that the type I peak persists even at 
a; = 16 meV with a frequency-independent peak width 
of 0.11 A~^. Based on this observation, we conclude 
that this mode is intrinsically inelastic at the intermedi- 
ate temperature regime (B). This conclusion is supported 
by the energy-scans with fixed Q at the type I AF Bragg 
point (Oil) in Fig. 2b. Furthermore, Q-scans at the type 
I point along [100] and [Oil] directions show that the 
magnetic scattering with > 10 meV around the type I 
Bragg point exhibits an isotropic peak in the (hOO)-(Okk) 
plane. We then observed a lower energy diffuse com- 
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ponent around the type II AF Bragg point, which has 
a much stronger intensity than the type I peak. The q 
width and scattering intensity of this component does not 
change with q along the hne connecting the (0,1.25,1.25) 
and (0.5,1.5,1.5) points within the limit of experimental 
errors. Therefore, we conclude that the magnetic scat- 
tering around the type II Bragg point is extended on the 
fee zone boundary. This component was concluded to be 
quasi-elastic by confirming energy scan at (1, -0.3, -0.3) 
which is close to the zone-boundary line. 
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FIG. 2. (a) Constant-energy scans along the [100] direc- 
tion at 32 K. Open and closed circles represent the data for 
hu) — 2 and 10 meV, respectively. The horizontal collimator 
is B-30'-60'-B with Ef = 14.7 meV. Absolute units for S{q,uj) 
were calibrated using the intensity of the transverse acoustic 
phonon. The right panel shows Reciprocal space for the fee 
lattice. The broken line shows the honeycomb zone boundary 
of the fee lattice in the (hOO)-(Okk) plane. Open circles, tri- 
angles and gray rectangles represent the AF type I, type II, 
and nuclear Bragg points, respectively. 

(b) Constant-Q scans at the type I AF Bragg point (Oil) and 
near the fee zone boundary (1,-0.3,-0.3) at 32 K. Open and 
closed circles show the data at (Oil) and (1,-0.3,-0.3), respec- 
tively. The data are divided by the square of the magnetic 
form factor at each position after subtraction of the back- 
ground. 

To extract information from the observed geometri- 
cal pattern of the critical scattering, we studied the 
momentum-transfer dependence of the energy integrated 



scattering law S{Q) = J dujS{Q,uj)/ \ F(Q) p ( F(Q) 
is the magnetic form factor), which is proportional to 
the real part of the static momentum dependent mag- 
netic susceptibility. Both above and below the Neel tem- 
perature T/vi the static susceptibility can be expressed 
in following form, x'(Q) - {(5 + (1 - /:(Q)/£(QAF)}~^ 
The parameter 6 specifies how much the actual temper- 
ature deviates from Tni, Qaf is an AF wave vector 
where the divergence of susceptibility at T = T/vi oc- 
curs, and C{Q) represents the Fourier transform of the 
interaction matrix. Considering the NN Ji, next NN 
J2 and next next (third) NN J3 AF interactions, one 
can express the Fourier transform of the interaction ma- 
trix£(Q) = - I Ji I {MQ) + R2MQ) + R3MQ)} in 
terms of the standard fee lattice geometrical factors (f>i 
and ratios R2 ~ J2/ Ji and i?,3 = J3 / Ji . 
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FIG. 3. Upper panel: Magnetic phase diagram of the fee 
lattice. Type I, type II and type III (Qaf = (1,1/2,0)) 
phases are indicated by I, II, and III, respectively. The area 
denoted by IC and the small triangle near the crossing of the 
lines separating I-II and II-III phases are regions where in- 
commensurate phases are found. Gray rectangles are possible 
locations of NiS2 in the phase diagram in the high symme- 
try phase. The arrow represents the possible change of the 
position in the diagram due to rhombohedral structural dis- 
tortion. 

Lower panel: Surface of critical scattering intensity for 
R2 = 0.434, R3 = 0.164, and 6 = 0.3. 

Since the value of the AF wave vector depends only on 
the parameters R2 and i?3, we analyzed the phase dia- 



3 



gram in terms of these coordinates (Fig. 3, upper panel). 
We found that the honeycomb structure (Fig. 3, lower 
panel) of the critical scattering in the (/lOO) — (Okk) plane 
is a general feature of the fee lattice with the AF interae- 
tions up to third NN. This pattern is robust in a rather 
large domain where i?2 < 0.6 and i?3 < 0.3. Therefore, 
to locate the position of NiS2 in the phase diagram, one 
have to study the fine features of the honeycomb struc- 
ture that are general for the fee lattice . To fix the 
values i?2 and i?3 for NiS2 in the high-symmetry phase, 
i.e. for T > T/v2, we calculated the critical scattering 
pattern at numerous points in phase I, and found three 
points (gray rectangles in upper panel of Fig. 3) that 
obey the following conditions obtained experimentally: 
(i) the scattering around the type I AF Bragg point is 
a symmetric circle-like peak; (ii) the contour plot of the 
scattering intensities around the type II AF Bragg points 
are ellipses aligned along the fee zone boundary; and (iii) 
energy integrated intensities at the type I and the type 
II AF Bragg points are equal within the limit of experi- 
mental errors. Therefore, one can conclude that NiS2 in 
a high-symmetry phase is situated in the vicinity of the 
gray rectangles and, hence, very close to the border sep- 
arating the I and II phases (see Fig. 3) This quasi- 
degeneracy also explains the anomalous temperature de- 
pendence of the spin fluctuations above Tni (regime A), 
namely the type I AF LRO temperature is suppressed by 
the frustrated magnetism, and the magnetic diffuse scat- 
tering is expected for a wide temperature range above 
Tni- 

Here, we discuss the fact that the two components of 
magnetic scattering have different properties in the in- 
termediate temperature regime (B). These components 
belong to different irreducible representations of the sym- 
metry group of the cubic lattice. 

Even when the type I AF LRO is established, both 
components are still present. Hence, the virtual ground 
state of the type II caused by magnetic frustration is 
quasidegenerate with the actual type I magnetic ground 
state of the cubic lattice. We emphasize that the 
quasidegeneracy occurs only accidentally between the 
two ground states which is incompatible in the cubic lat- 
tice symmetry (i.e. can be removed by a change of AF 
couplings) according to the Wigner theorem in group the- 
ory. Therefore, the system is not exactly frustrated but 
nearly frustrated. However, the nearly frustrated situa- 
tion can be resolved by the structural distortion which 
lowers the lattice symmetry. As a consequence, the type 
I and the type II AF LROs become compatible in the 
reduced symmetry group. In order to elucidate possi- 
ble consequences of SPT we studied small deviations of 
the magnetic exchange parameters constrained by the 
symmetry of the rhombohedral phase. Then we found 
that comparatively small changes of magnetic couplings 
(< 10~^) can induce a shift in the phase diagram denoted 
by the arrow in Fig. 4. A notable feature of the lowest 



transition at Tn2 is that the system avoids the phase of 
nearly frustrated magnetic interactions by lowering the 
lattice symmetry. 

In conclusion, in the paramagnetic and high symme- 
try phase (A) above Tni, NiS2 is a paramagnet where 
anomalies are governed by geometrical frustration and 
competing interactions. The second-order transition at 
Tni into the LRO type I AF state does not remove the 
nearly frustrated situation. The magnetic diffuse scatter- 
ing along the fee zone boundary, indicating the quaside- 
generacy of the type I and II AF LROs, persists down 
to the frustration-avoiding SPT at Tn2, where the type I 
and the type II phases probably become compatible due 
to the rhombohedral structural distortion. 
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